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The demixing of two-component fermions in optical lattices under a spin-dependent external 
potential is investigated using the density-matrix renormalization group method. The influence of 
on-site interactions (u) and the ratio of the spin-dependent external potentials (7) on the demixing is 
discussed. The j-u state diagram is numerically mapped out, from which we distinguish the phase- 
mixed region from the phase-separated one. A minimum point appears in the phase boundary. On 
the right part of the point an insulating core is formed in the spin-up densities during the demixing, 
while on the left one only metallic phase exists. 
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I. INTRODUCTION 



II. THE MODEL 



Cold-atom experiments provide a clean and versatile 
platform for simulating interacting quantum many-body 
systems and investigating new quantum phases p], [2j. 
In a system of cold atomic mixtures jsj, spin-dependent 
external potentials can be introduced by the different 
masses or the hyperfine states of the constituent Q. 
Spatial phase separation (i.e., demixing) of the compo- 
nents happens at repulsive interactions, which on the one 
hand will diminish the cooling efficiency in a sympathetic 
cooling procedure On the other hand, at stronger 
repulsive interactions and/or tighter confinement, this 
procedure provides a way to demix the different atomic 
species [6|-[sl- 

Different types of mixtures in optical lattices [9| have 
been studied, including boson mixture s II OIL fermion mix- 
tures of same or different species jlll-[l3|, which can be 
well modelled by a lattice model interacting through a 
Hubbard-type term with confining potentials. By detun- 
ing asymmetrically the laser frequencies with respect to 
the two hyperfine states |7j, or by trapping the two dif- 
ferent atoms of unequal masses, one realizes the different 
types of spin-dependent external potentials [l4]-[l6|. 

For a better understanding of the interplay of the spin- 
dependent external potentials and the repulsive interac- 
tion on the process of demixing, it is essential to have a 
complete state diagram, from which one can easily find 
the optimal parameters to realize the demixing or to con- 
trol the cooling efficiency. In this paper, we study the 
demixing of two fermion species (taking as pseudospins) 
of same masses in optical lattices with spin-dependent 
external potentials. 

This brief report is organized as follows. In Sec. II the 
modified Hubbard model under a spin-dependent exter- 
nal potential is introduced. Section III contains the main 
numerical results. Finally, we summarize the paper. 



A two-component Fermi gas in optical lattices with a 
spin-dependent external potential can be described by a 
modified one-dimension (ID) Fermi-Hubbard model as 



i,a i 

+ ^V a [i-(L-l)/2fn ia , 



(1) 
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where t is the hopping parameter between nearest neigh- 
bors and U the on-site repulsive interaction. The 
strength of the spin-dependent external potential is V a . 
The two components fermions are represented by pseu- 
dospin a =t, j. L is the length of the lattice, which is 
chosen to be long enough to keep the atomic densities at 
edges to be zero. c\ a (c^) is the creation (annihilation) 
operator on site i with spin a. The density operator for 
spin (7 is hifj. 

By taking t as the unit of energy, we rewrite the Hamil- 
tonian as 



H = ~XX 4x^+1^ +H.c.) ntfnu 

i,cr i=l 

+ ^^[i-(L-l)/2] 2 ( 7 n lt+ n li ) ) (2) 



with rescaled u = U/t and v a = V a /t. 7 = v^/v^ is the ratio 
between the spin- up and spin-down external potentials. 
For convenience, we assume that > v±. Thus, 7 > 1. 
The number of fermions of each species is N a = Y,i ^ia - 
Ei(ftio-) with riicr the spin-resolved density and the total 
number of fermions Nf = + N±. The polarization P = 
(N^-N\)/Nf. The local magnetization is defined as mi = 
(n^ -nn)/2. When 7 = 1, the model recovers the mostly 
studied Hubbard model in the presence of a harmonic 
trap [l7l-[2lj. We are interested in the effects induced 
by the spin-dependent external potentials. For atomic 
mixture of different species in optical lattices formed by 
lasers of wavelength 1024nm the ratio of the trapping 
potential 7 can reach 2.3 for fermionic mixtures of 6 Li 
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and 40 K, and 2.5 for 6 Li and 87 Sr. Larger 7 is hard to 
achieve at the present experiment but conceivable. 



III. NUMERICAL RESULTS 

We study in the following the effects of u and 7 on the 
phase separation. Our simulation is done by u sing the 
density-matrix renormalization group (DMRG) |22l l23j. 
In principle the number of states kept are 300 per block 
and 30 sweeps are used (when checking the numerical 
accuracy, up to 1000 states per block and 50 sweeps are 
used), with a truncation error in the final DMRG step 
usually smaller than 10~ 7 . Without mentioning in the 
text, we will choose the spin-unpolarized system with 
fixed particle number Nf = 40 and v± = 10~ 3 . 

We illustrate the physics induced by 7 in Fig. [TJ In 
panel (a), due to the same external potential, the spin- up 
and spin-down atoms have the same densities and accu- 
mulate at the potential minimum. With the increase of 
7, the spin- up atoms suffer from stronger trapping poten- 
tial, and part of the spin-down atoms are forced out of the 
system center, as the result of the competition between 
the potential energy and the interaction energy, which is 
shown in panel (b). Panel (b) including (a) is the case 
where the spin-up and spin-down atom mixture coexists 
in the center of the trap and forms a phase-mixed (PM) 
region. In panel (c), spin- up and spin-down atoms are 
becoming phase-separated at a critical j c = 5.1, decided 
when the local occupation for spin-down atoms disap- 
pears at the center (Numerically, we get the critical value 
of 7 C when ng & 10~ 3 with i ranging over the two cen- 
tral points j24j|.)- Beyond the critical j c , the spin- up and 
spin-down atomic clouds have diminished spatial over- 
lap and are phase-separated @, H| • A completely phase- 
separated case is shown at 7 » 7 C in panel (d). Panel (c) 
including (d) defines a phase-separated (PS) region. For 
comparison, in Fig. [TJ we present also the results based 
on the local spin-density approximation (LSDA) by us- 
ing a recent parametrized energy functional for the ID 
Hubbard model [25|. In the LSDA, the energy density of 
the inhomogeneous system is approximated at each spa- 
tial position by the one corresponding to a uniform gas 
at the local value of the density. We find that the LSDA 
gives qualitatively the same results as those of DMRG. 
However, the performance of the LSDA scheme at weaker 
interactions deteriorates with decreasing particle number 
where the hopping kinetic energy processes. As a result, 
the regions close to the edges of the trap and the phase- 
separation areas are also those where the LSDA is less 
accurate. 

To have a complete understanding on the effect of 7 
on the phase separation, in Fig. (2j the contour plot is 
shown for the spin-up density (up panel) and spin- 
down density n± (down panel) with 1 < 7 < 20. When 
7 > 12, an insulating core of fully polarized (spin- up) 
atoms (n^| = 1) is formed due to the strong confinement, 
which is incompressible and originates from the Pauli ex- 
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FIG. 1: (Color online) The density distribution of the two 
fermionic components under a spin-dependent harmonic po- 
tential v a (dashed line for v\ and dot-dashed for v\). The 
parameters are u = 5, N f = 40, and v± = 0.001. The DMRG 
data (open circle for and open square for n{) are compared 
against the LSDA ones (solid line for n\ and dotted line for 
n;). Panel (a): The two spin components have same distribu- 
tion due to the spin- independent external potential at 7 = 1. 
Panel (b): A part of spin-down atoms is forced out of system 
center with larger 7 = 1.4. Panels (a) and (b) belong to a PM 
region. In panel (c), at the critical point j c = 5.1 when there 
is no spin-down atoms left at the center of the system, a PS 
region starts to form. Panel (d): A real PS region is formed 
where the two components are demixed at 7 = 20. 




FIG. 2: (Color online) A contour plot for the influence of 7 on 
the demixing of the two components. The parameters are the 
same as in Fig.[T]but with continuously varying 7. Top panel: 
spin-up density (nt) and down panel: spin-down density (nt) 
are shown, respectively. With the increase of the 7, spin-up 
atoms become more accumulated and the spin-down atoms 
are forced out of the center of the system gradually due to 
the stronger repulsive interactions between them. 



elusion principle. As a result, the magnetization tends 
to 0.5. Around 7 ~ 20, stable density distributions are 
formed and no more changes can be observed by further 
increasing 7. 

The effects of the repulsive interaction u (with fixed 
7 = 3) on the phase separation are given in Fig. [3l For 
a non- interacting system (u = 0), the two components 
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FIG. 4: (Color online) (Left panel) Ground-state energy Eqs 
(open circle) as a function of 7. The parameters are the same 
as Fig. El The energy without the contribution from the spin- 
down external potential Ea (open square) is also plotted. The 
dash line is the linear fit for Eqs with Eqs = -24 + 0.737- 
(Right panel) Ground-state energy Eqs (open circle) as a 
function of u. The parameters are the same as Fig. [3] The in- 
teraction energy E u is also shown (open square) , which starts 
from zero when there is no interaction, reaches a maximum 
at U = 5, and goes to a constant when the two components 
segregate almost completely. The dash line is a constant at 
E = -28.22. 



FIG. 3: (Color online) The influence of u on the demixing 
of the two components. The system consists of iVf = 20 and 
iVj = 20 fermions. The spin-down trap strength is 10~ 3 and 
7 = 3. From top to bottom, the total density distribution 
(n), spin- up density (rif), spin-down density (ri|), and the 
magnetization (m) are shown, respectively. 



keep mixed. The spin-down atoms are forced out of the 
center due to the strong repulsive interaction between the 
components when increasing u. The system reaches the 
critical PS point at u c = 8. The two fermion components 
are fully separated around u ~ 15 at which the density 
distributions become stable. 

The interplay of 7 and u segregates the two compo- 
nents. The spin- up atomic cloud that suffers from a 
stronger external potential tends to accumulate at the 
center. As a result, the spin-down atoms are repelled 
from the center due to a larger repulsion, which decreases 
the interaction energy but increases the external poten- 
tial energy. The equilibrium density distribution is the 
result of the balance between these two opposite effects. 
The shape of the density profiles becomes stable at fur- 
ther increasing 7 and u. To understand this point, we 
plot in Fig. [4] the energies as functions of 7 and u, re- 
spectively. In the left panel, we study the ground-state 
energy 

E GS = T + uD + E l ext+1 El xt , (3) 

and the energy except the last term Ea - T + uD + 
E^ xt as a function of the ratio 7. Here T is the ki- 
netic energy, D the double occupancy D = ££i(^Wt^l)> 
and the respective external potential energy is E° xt = 

v i - ( L ~ (^>- We notice that at 7 > 20, 



Egs becomes a linear function of 7: Eqs - -24 + 0.737 
and Ea goes to a constant -24, which indicates that the 
densities are becoming stable and the only contribution 
to the energy comes from the last term linear to 7. In the 
right panel, we study the ground-state energy Eqs and 
the interaction energy E u - uD as a function of the in- 
teraction strength u. We can see that at u > 22, Eqs and 
E u become constant, which indicates that the two com- 
ponents are demixed almost completely except the tails 
of the clouds. Under the stable density distributions, the 
on-site interaction contributes a constant energy. 

To have a overview on the transformation of the sys- 
tem from PM to PS region under the interplay of the ra- 
tio of the external potential and the repulsive interaction, 
we show the 7 - u state diagram in Fig. [5] for different 
parameter systems. For the unpolarized systems, three 
different v i = 0.001,0.003, and 0.006 are shown. A po- 
larized system of v± = 0.001 and P = 0.5 is also included. 
As mentioned in the Ref. [H , there exists a critical repul- 
sive interaction strength (u c ) above which the demixing 
evolves. All the critical points form the phase boundary 
which distinguishes the PM phase from the PS one. For 
the unpolarized system, the different v± give very simi- 
lar phase boundaries. For the systems of the same ratio 
7, the demixing of the polarized system is easier due to 
more spin- up atoms and less spin-down atoms in the trap. 
Thus, smaller u c is observed for the onset of the phase 
separation. 

By carefully investigating the different phase bound- 
aries, we find that there is a minimum point in each of 
the phase boundaries, labeled by (70, uo), which is closely 
related to whether the insulating core of spin-up atoms 
is formed in the trap center. To understand this mini- 
mum point, we show in Fig.[6]the 3D plot of spin-up and 
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FIG. 5: (Color online) The 7 - u state diagram for different 
systems. For the unpolarized systems, the phase boundaries 
for three different v± = 0.001,0.003, and 0.006 are shown by 
down- triangle, circle, and square, respectively. For the po- 
larized system for v± = 0.001 (P = 0.5), it is shown by the 
up-triangle. The number 2 and 3 are the corresponding min- 
imum position in the phase boundary, while 1 and 4 mark 
only the corresponding minimum positions in the w-axis. The 
dashes lines are just a guide for eyes. 
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spin-down densities for the system of Nf = 40, P = 0, and 
v± = 0.006. 7 is fixed to be 3 and 5 in the upper and bot- 
tom panels, corresponding to the left and right part of the 
minimum point [i.e., the third point: 3 (70 « 4,ii « 4.32) 
in Fig. [5] in the phase boundary, respectively. We check 
the spin-up and spin-down densities by varying the inter- 
action strength. We find that the local spin-up densities 
for 7 = 3 are always smaller than one (n^ < 1) during 
the demixing process, i.e., there is no insulating plateau. 
While the spin-up densities for 7 = 5 form an insulat- 
ing core (n^i = 1) during the demixing process. From 
the above analysis we know that, the larger the 7, i.e., 
the tighter the confinement, the easier is to deplete the 
spin-down atoms. As a result, smaller u c is needed to 
achieve PS. This is the reason why the phase bound- 
ary firstly goes down when increasing 7. However, when 
the insulating plateau appears, the depletion of the spin- 
down atoms becomes more difficult and stronger on-site 
interaction is needed to achieve PS. That is, larger u c is 
needed, which makes the phase boundary go up slightly. 



IV. CONCLUSIONS 

In this paper, we discuss the demixing of fermion mix- 
tures in optical lattices under the presence of a spin- 
dependent external potential by means of the DMRG and 
LSDA methods. The effects of the ratio of the spin-up 
and spin-down dependent external potential 7 and the 
repulsive interaction u on the phases of the system are 



FIG. 6: (Color online) The 3D plot for the spin-up and spin- 
down densities as functions of the site i and the interaction 
strength u. The parameters chosen are, Nf = 40, P = 0, 
and vi = 0.006. The upper panels is for fixed 7 = 3, locating 
at the left of the minimum point in the phase boundary with 
u e [4, 9]. In this case, the spin- up densities are always smaller 
than one during the demixing process. The Bottom panels is 
for fixed 7 = 5, locating at the right of the minimum point in 
the phase boundary with u e [3, 4.6]. In this case, the spin-up 
densities form an insulating core during the demixing process, 
which leads to an extra energy needed to deplete completely 
the spin-down atoms. 



analyzed in details. The competition between the spin- 
dependent potential energy and the repulsive interaction 
energy is responsible for the phase separation. The ex- 
ternal potential tends to push the atoms to the center of 
the system, while u tends to delocalize atoms. The j-u 
state diagram is obtained in which the phase boundary 
distinguishes the PM from the PS phases. A minimum 
point in each of the phase boundaries is observed. In 
the left part of this point, the local spin- up densities 
are always smaller than one which makes the depletion 
of the spin- up comparatively easier. While in the right 
part of this point, the spin- up densities form an insu- 
lating core (riij = 1) during the demixing process. The 
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insulating plateau makes the depletion of the spin-down V. ACKNOWLEDGEMENTS 

atoms more difficult and as a result, the phase boundary 
in this case goes up slightly. 

The detailed state diagram found in this paper is help- 
ful in choosing the optimal parameters to achieve demix- 
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